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We study the effect of the transverse degrees of freedom on the velocity of 
sound in a Bose-Einstein condensate immersed in a one- dimensional optical 
lattice and radially confined by a harmonic trap. We compare the results of 
full three-dimensional calculations with those of an effective ID model based 
on the equation of state of the condensate. The perfect agreement between 
the two approaches is demonstrated for several optical lattice depths and 
throughout the full crossover from the ID mean-field to the Thomas Fermi 
regime in the radial direction. 



Ever since the achievement of Bose-Einstein condensation (BEC) with 
atomic vapors, an extensive activity in the study of the excitations of these 
systems has been carried out. A particular example is the exploration of 
sound excitations in elongated condensates l ^ l '^ l '^ l ^ l ^ l The propagation of 
sound in one-dimensional (ID) optical la ttices has b een the subject of re- 
cent studies performed both in the linear I I I I I I I and nonlinear regimes 
To capture the main physics, the system can be conveniently described 
by means of effective ID theories, which account for the transverse degrees 
of freedom through a renormalization of the coupling constant. Still, it is 
important to check the effect of the transverse degrees of freedom since they 
are known to play a relevant role in some cases 

In this paper, we show that the sound velocity in presence of a ID 
lattice and transverse harmonic confinement can be obtained by employing 
an effective ID approach, and demonstrate that the effect of the transverse 
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harmonic confinement is completely accounted for by the expression 



dN 



where is the number of atoms per lattice well and m* is the effective 
mass. The effect of lattice and transverse confinement both on the equation 
of state ^(A^) and m* is discussed in detail. 

Let us consider an axially symmetric condensate, which is radially con- 
fined by the potential T^io = rnuj'j_r'^ /2, and immersed in a ID periodic po- 
tential of the form Mat (■2) = s Eji sin^ (vrz/d), where d is the lattice period, 
Eji = q^/2m the recoil energy, qs = h-w/d the Bragg-momentum, and s the 
lattice depth in units of the recoil energy. We consider any additional ax- 
ial confinement to be negligible, so that the Gross-Pitaevskii (GP) equation 
takes the form 



2m 



^'(r, z) = /i^(r, z) , 



where is the number of atoms per lattice well and the condensate wave- 
function is normalized according to J d'^r J^^^^^dzl^lr, z)\'^ = 1. 

In the linear regime, the small oscillations around the condensate ground 
state can be described by expanding the order parameter as ^{r,z,t) = 
g-iMt/?i[^(r^ z) + J2juq Uj^ge-''^^-^'}'^^ + Wj%e''^^''(9)*], where the functions uj^g 
and Vjuq are the Bogoliubov quasiparticle amplitudes of quasimomentum q, 
and quantum numbers j, v representing the band index and the number of 
radial nodes respectively. In particular, the lowest energy solutions, with 
J = 1 and = 0, are axial phonons that propagate with the sound velocity 
c. The latter is defined as the slope of the lowest branch in the excitation 
spectrum Uji,{q), through the relation c = hdLOio{q) /dq\q=o. Therefore it 
can be obtained by numerically solving the spectrum of the full 3D system, 
as discussed in^^. 

In the following we will show that the system can be equivalently de- 
scribed by an effective ID approach. The basic assumption we make is that 
the condensate wavefunction can be factorised into a longitudinal and a ra- 
dial part ^(r,z) = ip{r)ip{z). The physics underlying this assumption is 
that interactions mainly affect the radial wavefunction c/j(r), while the axial 
wavefunction tp{z) is that of a single particle in the periodic potential. This 
factorization is exact for s = 0, where the system is uniform in the longi- 
tudinal direction. In the case s 7^ 0, the solution of the GP equation for a 
transversally uniform system shows that for the densities of many current 
experiments^^^E^^the equation of state and the effective mass can be cal- 
culated neglecting the density dependence of ip{z)^. The remaining relevant 
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density dependence of these quantities arises from their exphcit dependence 
on the number of particles (see below). The range of densities for which 
the factorization ansatz is correct becomes larger as the lattice depth s is 
increased. 

Under the factorization assumption, and using the normalisation con- 
dition for the wavefunction, one can write two effective GP equations 



2m 




z)\'^dz N\(p{r)f 



V9(r)|Vr N\ij{z)\ 



i,{z)=,,Mz), (4) 



where is the number of particles in each well at equilibrium, fi^ = ^ho 
and fi^ = n - eiat, being eho = J d'^rip*{r) -^Vf. + Vho{r) (fir) and eiat 
the single particle ground state energy of the periodic potential. 

The effective mass is calculated from the lowest energy Bloch band 



e{k) = / rk{z) 

J-d/2 



|;v?+M..(.)+|^fe(.)r^ 



^kiz)dz, (5) 



according to the general relation 1/m* = d'^e/dk'^\k=o In Eq.© the 
explicit dependence of the energy per particle on density appears through 
the renormalized number of atoms N = Na^ J \ip{r)\^d'^r, where a± = 
^/h/mu>±. Consistently with the above discussion, the Bloch state ipkiz) 
at quasi-momentum k is taken as the single particle solution of @. Then 
Ipkiz) is independent of the radial confinement and the effect of the radial 
trapping on m* is captured by its dependence on the renormalized number 
of atoms A^. 

Equation (jSJ describes a EEC uniform along z and harmonically con- 
fined along r (infinite cylinder). The effect of the lattice is included in the 
effective coupling constant 

/■d/2 

g = gd mz)\''dz, (6) 

J-d/2 

where again ip{z) is the single particle ground state of @. Without radial 
confinement, in a system with 3D average density n, the equation of state 
IJ-±{n) = gn is linear in the density, since g depends strongly on the optical 
lattice depth, but not on density. As mentioned above, this result is known 
to be correct in the presence of a lattice for sufficiently small densities ^. 

In order to discuss the effect of the transverse degrees of freedom, we 
concentrate now on the radial dynamics. The presence of the lattice is de- 
scribed by the effective coupling constant ^, so that the crossover from frozen 
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radial dynamics (ID mean-field) to the Thomas-Fermi (TF) regime, can be 
described in full analogy to the case of the infinite cylinder^. Introducing 
the dimensionless variable p = r/a± and defining f{p) = a±ip{r), Eq.® can 
be rewritten in the dimensionless form 

where uid = N/d is the linear average density, a is the effective scatter- 
ing length obtained from the relation g = ATrh^a/m, and / is normalized 
according to / d?p\f\'^ = 1. Notice that the solution of the above equation 
depends only on the parameter dniu, yielding an exact scaling behavior for 
the radial wavefunction as a function of the transverse trapping. 

The dependence of the chemical potential p±/hu;± on dniD determines 
the equation of state in the presence of the radial trap. For dniu ^ 1, 
the system is in the TF regime and one finds fi±/hujj_ = 2{dniD)^^'^ ■ In 
the opposite limit dniD <C 1, instead, the radial profile is gaussian and 
the equation of state is given by pL_\_/hoj±_ = 1 + 2dniD ■ Given the solution 
p±/huj± as a function of the universal parameter dnnj, the velocity of sound 
can be calculated according to Eq.(^. This procedure yields the result 
c = y^gn{0)/2m* both in the TF regime {dnm ^ 1) and in the ID mean- 
field limit {dniD 1) To find these limiting results, we have used the 
fact that niD = {4:n'^h'^d/m^io'j_)n^{0) in the TF regime and nm = 7Ta']_n{0) 
in the ID mean-field limit. 

The result for the sound velocity in the intermediate regimes can be cal- 
culated numerically, just based on the knowledge of the following quantities: 
(i) the renormalized scattering length d for any given value of s, obtained 
from a ID calculation of the single particle ground state solution in presence 
of the lattice; (ii) the equation of state and for the infinite cylinder at any 
dniD', (iii) the effective mass m* at any given value of s for the renormalised 
number of atoms per well A^, obtained from Eq. © by a ID calculation of 
the single particle Bloch solutions. This procedure leads to the results shown 
in Fig^ where the theoretical predictions are compared with the results of 
the full 3D numerical calculations for several values of the lattice depth. The 
agreement is perfect in the full crossover. 

Furthermore, it is interesting to note that by relating nid to the density 
at the center n(0) one can show that the value c = y^gn(0)/2m* provides an 
estimate for the sound velocity which is correct within 5% throughout the 
whole crossover from TF to ID. 

We also notice that in the TF regime, driiu ^ 1, the sound velocity 
in presence of transverse harmonic confinement corresponds to that of a 
transversally uniform condensate with the average density n, c = y^gn/2m* . 
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Fig. 1. Velocity of sound as a function of dniD, for three values of the 
lattice depth (from top to bottom s = 0, 5, 10). Results of full 3D numerical 
calculations (points) are compared with the prediction derived from the ID 
equation of state (lines). Note that d/a = 1.49 and d/a = 1.98 for s = 5 and 
s = 10 respectively. 

In fact, in this regime the radial profile is an inverted parabola, leading to an 
aver age d ensity n = J dxdyn{x,y) /irdR^ equal to half of the central density 
j^(0) HIM], This explanation in terms of n works well in the TF regime, 
and to some extent also in the ID mean-field regime, where the definition 
for the average density is replaced by n = / dxdyn?{x^y) / J dxdyn{x,y)^. 
However, since the concept of average density is not well defined in general, 
it does not provide a general criterion for the determination of the sound 
velocity in the intermediate regimes. 

In conclusion, we have demonstrated that the equation of state of an 
infinite cylinder captures the effect of the transverse degrees of freedom on 
the sound velocity in a condensate immersed in a ID optical lattice, upon 
a suitable renormalization of the mass and of the coupling constant. This 
result is obtained by means of a factorization ansatz for the condensate wave- 
function. The predictions of the factorized model yield perfect agreement 
with full 3D numerical calculations in a wide range of optical lattice depths 
and in the full crossover from the ID mean field regime to the TF limit. 
The factorization ansatz rely on the assumption that in absence of radial 
trapping the equation of state is linear in the number of particles. This 
condition is always fulfilled at s = and in the tight binding regime, while 
for intermediate s it is satisfied to a good approximation at experimentally 
relevant densities. In any case, Eq.(p3) is completely general, and could be 
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easily checked also beyond the limit of validity of the factorization ansatz. 
In that case the determination of n{N) would require from the beginning a 
full 3D calculation. Instead, when the factorization assumption holds, the 
simple model presented in this paper can be used to extract the correct value 
of the sound velocity. 
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